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Abstract. We report the detection of 1 x 10 7 Mq of molecular gas in the central region of the S0/E7 galaxy 
NGC4550, inferred from observations of CO(l-O) emission. Dust is detected in HST WFPC2 images and found 
to be asymmetrically distributed around the nucleus, only extending to a galactocentric distance of 7" (600 pc). 
The shape of the CO emission profile is consistent with a molecular gas distribution following the dust. The 
distribution of the dust and gas in the center could be the result of an m — 1 instability, which is the fastest 
growing unstable mode in counterrotating stellar disks. On a global scale the molecular gas in NGC 4550 is stable 
against gravitational collapse but nevertheless star formation appears to be ongoing with normal star formation 
efficiency and gas consumption time scales. The stellar velocity dispersion in NGC 4550 resembles that of elliptical 
galaxies. It is therefore likely that a hot X-ray emitting plasma limits the lifetime of the molecular gas, that must 
arise from a recent (<Cl Gyr) accretion event. 
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1. Introduction 

While counterrotating stellar systems are quite common 
among elliptical galaxies, counterrotating stellar disks are 
very rare. Despite dedicated searches (cf. Kuijken et al. 
1996), there are only two known galaxies with large coun- 
terrotating stellar disks; NGC 4550 (Rubin et al. 1992) 
and NGC 7217 (Merrifield & Kuijken 1994). There are a 
few cases which do contain counterrotating stellar disk 
components, but where the secondary component is much 
less extended than the primary one and is confined to the 
central region: NGC 3593 (Bertola et al. 1996), NGC 4138 
(Jore et al. 1996) and NGC 7331 (Prada et al. 1996). On 
the other hand Bertola et al. (1992) and Kuijken et al. 
(1996) report that ~20-25% of all gas disks found in SO 
galaxies counterrotate with respect to the stars. If this gas 
would form stars in such numbers, that more than 5% of 
the stellar disks were counterrotating, a much higher frac- 
tion of galaxies with counterrotating stellar components 
would be observed. 

The most prominent example of a galaxy with two 
counterrotating stellar disks is NGC 4550, an S0/E7 
galaxy in the Virgo cluster, which contains two expo- 
nential stellar disks with approximatively the same scale 
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length and central surface brightness (Rix et al. 1992). 
One of the disks also contains a small amount of ionized 
gas (Rubin et al. 1992). The stellar disk corotating with 
the ionized gas has a slightly higher velocity dispersion 
and lower maximum rotational velocity than the stellar 
disk counterrotating with respect to the gas (Rix et al. 
1992). Recent results by Faundez et al. (2000) also suggest 
that the disk containing gas is sligthly more massive, more 
vertically extended and dynamically hotter than the other. 
The equivalent widths of the Ca H and K features are the 
same for the two disks and similar to those found in ellip- 
tical galaxies, suggesting that the stellar disks are old and 
approximately coeval (Rix et al. 1992). NGC 4550 con- 
tains a low luminosity active galactic nucleus (AGN), seen 
through broad (^10 4 kms _1 ) and double-peaked emission 
lines (Ho et al. 2000). 

Is the origin of the massive counterrotating stellar 
disk in NGC 4550 different from that of the numerous 
less massive counterrotating gaseous (and in some cases 
stellar) disks seen in other galaxies? Misaligned angular 
momentum axes are considered to be clear signatures of 
merging (cf. Thakar & Ryden 1996; Thakar et al. 1997). 
Strong disk mergers usually lead to the destruction of the 
disks and the formation of elliptical systems. Recently, 
however, Pfenniger (1999) has shown that similarly sized 
disks, but with oppositely aligned angular momentum vec- 
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tors may merge without significant heating of the stellar 
populations. Here we report the detection of molecular 
gas and present optical images of the dust distribution 
in NGC 4550. Consequences for the origin of this gas and 
dust are evaluated. 



2. Observations 

2.1. Molecular data 

The CO observations were made with the IRAM (Institut 
de Radio Astronomie Millimetrique) 30-m telescope on 
Pico Veleta in Spain in October 1992 and May 2000. The 
full width to half power beam size was 9b = 23" (1.9 kpc 
at a distance of D = 16.8 Mpc; Tully 1988). 3mm SiS 
receivers were used to observe the CO (1-0) line with a 
nutating subreflector, switching symmetrically ±240" in 
azimuth with a frequency of 0.5 Hz. Typical system tem- 
peratures were 400 K (Oct. 1992) and 350 K (May 2000) 
on a main beam brightness temperature (T m b) scale. In 
1992, filterbanks were used with a channel spacing of 2.6 
kms -1 and a total bandwidth of 500 MHz. In 2000, a dig- 
ital spectrometer was employed; the channel spacing was 
3.2 kms -1 and the bandwidth was 640 MHz. Final spectra 
were obtained by adding all individual spectra weighted 
with the inverse of the square of their rms noise, remov- 
ing a first order baseline and binning the channels to an 
effective velocity resolution of 20-31 kms -1 . 

Typical pointing corrections determined by regular ob- 
servations of nearby continuum sources were 3-4". The 
temperature scale used here is T m b = T£ x F c g/B c s, with 
the forward hemisphere efficiency F e s — 0.95 and the main 
beam efficiency B c s = 0.76. In 1992, image sideband rejec- 
tions were assumed to be 7dB, but were not determined 
for each individual pointing. An error in the sideband re- 
jection of only ±ldB would cause a calibration error of 
±5%. In 2000, image sideband rejections were ~25dB; 
therefore uncertainties in the image sideband rejection do 
not significantly contribute to the total calibration error. 
As a consequence, the temperature scale of these measure- 
ments is more reliable. 

2.2. Optical data 

We retrieved Hubble Space Telescope (HST) archival data 
of NGC 4550 obtained during Cycle 4. The data consist of 
Wide Field/Planetary Camera 2 (WFPC2) images with 
the F555W (~V-band) and F814W (~I-band) filters. The 
data were obtained on December 10, 1994, (proposal 5375, 
PI: V.C. Rubin), with exposure times of 3x400 seconds for 
each filter. The WFPC2 data were processed through the 
standard pipeline where bias, dark, and flatfielding correc- 
tions were performed and photometry keywords were cal- 
culated. Subsequent processing was done using standard 
tasks in IRAF/STSDAS. Multiple images were simulta- 
neously co-added and cosmic-rays removed. The images 
shown here are from the Planetary Camera (PC) with 



a pixel scale of 0.046", corresponding to 3.8 pc at D = 
16.8 Mpc. 



3. Results 

3.1. Molecular gas 

In Fig. 0a we show the CO(1-0) emission profile observed 
in May 2000 towards the center of NGC 4550. A weak but 
clear signal is detected. The emission extends from —225 
to —535 kms -1 , which is similar to that of stars and the 
ionized gas (Rubin et al. 1992; Rix et al. 1992). The veloc- 
ity integrated intensity is 0.9±0.1Kkms _1 . With a A^h 2 - 
to-2co conversion factor of 2.3 x 10 20 cm~ 2 (Kkms -1 ) -1 , 
and D = 16.8 Mpc (Tully 1988), this corresponds to an H 2 
mass of (1.3±0.1) x 10 7 M Q . The total stellar mass can be 
estimated from the i?-magnitude, assuming an M/L ratio 
of 5, to be 2 x 10 10 M . This gives a gas mass fraction of 
only -0.1%. 

The CO emission profile is skewed relative to the 
systemic velocity of the galaxy (V sys ~ 380 kms , e.g. 
Huchtmeier & Richter 1989; Rubin et al. 1997), being 
stronger at higher velocities. Dividing the profile into a low 
and high velocity part, measuring the velocity integrated 
intensities between 220-380 kms -1 and 380-520 kms -1 
we get {Iugh - how/Irms) ~ 1.9. For a molecular gas dis- 
tribution being symmetric with respect to the minor axis 
crossing the nucleus, we expect a symmetric emission pro- 
file. Skewed profiles like the one seen in Fig. [j]a can then be 
caused by a pointing offset of the telescope. However, in 
this case we observed the same position for several hours, 
with pointing checks at regular intervals. Small pointing 
uncertainites should therefore result in a general broaden- 
ing of the telescope beam resulting in a symmetric emis- 
sion profile rather than a skewed one. Also, during the 
course of the observations the source rotated with respect 
to the telescope beam counteracting any systematic point- 
ing offset. It is therefore likely that the skewed shape of 
the emission profile in Fig. [j]a is real and reflects an asym- 
metric distribution of the molecular gas. 

In 1992 we observed CO(1-0) emission towards 6 po- 
sitions in the central region of NGC 4550. Offsets were 
(Aa, AS) = (0",0"), (0",+10"), (0",+20"), (0",-10"), 
(+10",0") and (-10",0"), relative to the center position 
(see Fig. IT]). No emission was directly detected above the 
3<7 level in any individual position but by convolving the 
observed data with a telescope beam of 0b = 36", we 
detect a clear signal. Omitting the northernmost posi- 
tion from the convolution we get the CO profile shown 
in Fig. |l|b. The exclusion of position (0",+20") does not 
change the overall appearance significantly. The emission 
profile is similar to the one shown in Fig. [l]a, i.e. the pro- 
file is skewed. The integrated intensity is higher in Fig Jib, 
but in view of the calibration uncertainties (see Sect. |2.l|) 
and the small number of measured positions that could be 
used for the convolution, the temperature scale in Fig.[l]a 
is more reliable. 
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Fig. 1. a) The CO(l-O) emission profile from NGC 4550 obtained with a full width to half power beam size of 6*b 
= 23" in May 2000. The velocity resolution is 20kms _1 . b) The CO(1-0) emission obtained by convolving the 1992 
data with a beamsize Oh of = 36" (see text). The velocity resolution is 31kms _1 . In both cases the telescope beam is 
centered on the nucleus of NGC 4550, at a = 12 h 35 m 30.9 s , 5 = +12° 13' 17" (J2000). The velocity scale is heliocentric 
(Vlsr. = Vhel + 4.3 km s -1 ) and the intensity is given in units of main beam brightness temperature. 



In order to test how sensitive the shape of the emission 
profile is to offsets of the gas from the nucleus, we made a 
simple model of a highly inclined (i — 75°; cf. Tully 1988) 
molecular gas distribution with the major axis along a 
north-south direction (for NGC 4550, PA = 358°; Rubin 
et al. 1992) and the receding part of the galaxy to the 
north (Rubin et al. 1992). For the molecular gas distribu- 
tion we assumed a gaussian distribution with a full width 
to half power source size of 9 S = 14" (see Sect. 3.3). For 



the velocity field we adopted a parametrization of the ro- 
tation curve resembling the one presented by Rubin et al. 
(1992). The resulting iso- velocity contours, 'observed' with 
a telescope beam size identical to the one used in the real 
observations (9b = 23"), are shown in Fig.0. The exact 
shape of the rotation curve is not crucial for demonstrat- 
ing the effect from offsetting the center of emission. In the 
top panels we show the case where the gas is centered on 
the nucleus. The observed spectra are symmetric, with a 
top-hat shape. This results from the fact that the gas is 
sampled mainly on the rising part of the rotation curve. 
In the lower panels the gas distribution has been shifted 
3" north, while keeping all other parameters fixed. The re- 
sulting emission profile is strongly skewed towards higher 
velocities and resembles the observed CO(1-0) profiles in 
NGC 4550 (Fig. |l]). 

Centering the convolution of the CO profiles observed 
in 1992 5" north and south of the nuclear position, the 
northern profile shows an integrated intensity that is 
30%±10% larger than the corresponding profile from the 
southern position. While this result is only marginally sig- 
nificant, it is nevertheless consistent with a skewed spatial 



distribution; most of the CO emission arises from locations 
north of the center. 



3.2. Atomic gas 

There are conflicting reports in the literature regarding 
21cm HI observations. While a few early observations re- 
ported detections with integrated fluxes between 1-9 
Jykms^ 1 others, including all recent observations, have 
resulted in non-detections. Using the Arecibo antenna 
DuPrie & Schneider (1996) report a non-detection at an 
rms level of 1.5 mJy. At a distance of 16.8 Mpc this corre- 
sponds to a 3ct upper limit to the HI mass of 7 x 10 7 M Q . 
Although this limit is considerably higher than the de- 
tected molecular gas mass, it shows that the total amount 
of gas in NGC 4550 is < 10 s M Q . 



3.3. Dust 

The asymmetric distribution of the molecular gas indi- 
rectly inferred by the lineshapes of our CO profiles (Fig.|lJ) 
can be traced directly by the dust distribution. Dust is 
visible as dark patches in the V-band image and becomes 
prominent in the V— I image (Fig. The dust can be 
followed to a radius of ~7" and is distributed in several 
arclets, reminiscent of spiral arms. The main morphologi- 
cal feature is, however, the non-axisymmetric distribution, 
with the northern part being significantly more prominent 
than the southern part. This could in principle be caused 
by a bar-shaped dust distribution seen at an angle to the 
plane of the sky, with the northern side being closer to 
the observer. The difference in background light behind 
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Fig. 2. Illustration of the effect of shifting the center of a molecular gas distribution relative to the galactic nucleus. 
In this model we 'observe' a gaussian shaped CO distribution with a full width to half power of 9 S = 14" with a beam 
size of Ob = 23". The rotation curve is taken from Rubin et al. (1992) and the inclination is set to 75°. The left panels 
show the CO distribution in grey scale and the iso- velocity contours in steps of 20kms _1 . The systemic velocity is 
set to 380 km s . The right panels show the resulting CO emission profiles. In the top panels we show the result by 
centering the CO distribution on the nucleus. In the bottom panels we shift the CO distribution 3" to the north, while 
keeping all other parameters constant. The resulting CO emission profile shows a distinct skewness. 



the dust would yield a higher contrast (see Fig. 3c) in the 
north. The extreme difference of the dust obscuration be- 
tween northern and southern side of the nucleus appears, 
however, to be too large to be explained in this way. An 
intrinsically non-axisymmetric distribution is more likely. 

NGC 4550 was detected by IRAS in the 60 and 100/zm 
bands (140 ± 31 and 220 ± 80 mJy). With D = 16.8 Mpc 
the total far-infrared (FIR) luminosity amounts to (1.0 ± 
0.3) x 10 8 L Q . The FIR flux corresponds to a dust mass 

M dust = 4.8xl0- 11 ^^- M , (1) 

where B v is the Planck function evaluated for a dust tem- 
perature Td = 39 ± 7K (60/100yum color temperature), 
Dm P c is the distance measured in Mpc (16.8 Mpc), S v is 
the flux measured in Jansky (here we use the 100/un flux 
of 220 ± 88 mJy) and k v is the mass opacity coefficient 
for which we used 2.5 m 2 kg -1 (Hildebrand 1983) and as- 
sumed k v oc u +1 . The dust mass is (1.2 ± 0.9) x 10 5 M©. 
The large uncertainty is mainly due to the uncertainty 
in the dust temperature. The corresponding gas-to-dust 



mass ratio is ~ 150±114, when incorporating a primordial 
abundance of helium. The ionized and atomic gas masses 
are assumed to be negligible in comparison to the molec- 
ular gas mass (cf. Rubin et al. 1992; DuPrie & Schneider 
1996). Thronson & Telesco (1986) found a typical gas-to- 
dust mass ratio of ^700 for normal spiral galaxies, which 
reflects the fact that the IRAS bands are not sensitive to a 
dust component colder than ^20 K. The lower gas-to-dust 
ratio found in NGC 4550 indicates that a relatively warm 
dust component dominates the far infrared continuum. 

The heating source of the dust can be (i) the ambient 
stellar radiation field, (ii) massive young stars embedded 
in the dust clouds, (iii) photons from a hot X-ray emitting 
plasma or (iv) collisions with electrons from a hot X-ray 
gas component. 

The observed visual magnitude within a 30" aper- 
ture centered on the nucleus of the galaxy is V=12.31 
(Sandage & Visvanathan 1978). This corresponds to an 
average surface brightness of 19.4 magarcsec~ 2 . Using a 
solar absolute magnitude My Q = 4.83 mag we get an av- 
erage intensity for the ambient stellar radiation field of 
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/* = 1.1 x 10 -2 ergs -1 cm~ 2 stcr -1 . The heating rate 
per dust grain by absorption of ambient stellar light is 
then Airl* Qabs ^a 2 « 1.4 x 1CP 11 ergs -1 , where we have 
used an absorption coefficient of 0.36 and a typical dust 
grain radius of 0.1//m (cf. Jones & Merrill 1976; de Jong 
et al. 1990). The equilibirum dust grain temperature can 
be estimated by equating the heating and cooling rates 
and becomes ~25K. This is 2a lower than the observed 
39 ± 7 K and suggests an additional heating source for the 
dust. Although there are no obvious signs for massive stars 
in NGC 4550 (cf. Fig. ||), the gas mass derived from the 
CO emission and from the FIR luminosity are consistent 
with each other assuming a normal gas-to-dust mass ra- 
tio. Hence it cannot be excluded that youn g ma ssive stars 
are embedded in the molecular gas. In Sect. [O] we will see 
that X-ray heating, either through X-ray photons or hot 
electrons, is inefficient in comparison with heating by the 
ambient stellar radiation field. 

3.4. X-ray properties 

NGC 4550 is undetected at X-rays. Einstein data give 
L x < 4.1 x 10 40 ergss -1 (Fabbiano et al. 1992). It is never- 
theless likely that a hot X-ray emitting plasma is present 
at some level. Stellar mass loss from evolved stars will 
contribute ~ 0.15 (L B /10 10 L s ) M yr _1 of gas to the in- 
terstellar medium (cf. Faber & Gallagher 1976; Sarazin 
1990). The large stellar velocity dispersion in NGC 4550, 
caused by the counterrotation, will heat this gas to X-ray 
temperatures. The gas lost from evolved stars will have an 
energy per unit mass of 3fcT/2/im p , similar to the energy 
per unit mass of stars, which is of the order <r 2 /2 (using 
the one-dimensional velocity dispersion since stellar orbits 
are confined to a plane). After thermalization this leads 
to an equilibrium temperature of (cf. Sarazin 1990) 



(im p <j 
3k 



2x 10 6 f T Y K 

VaOOkms- 1 / 



(2) 



An additional source of heating comes from supernovae, 
but due to the high relative velocity dispersion cr* in 
NGC 4550, it is likely that energy input from stellar mass 
loss will dominate the heating. Assuming that the gas re- 
mains gravitationally bound to the galaxy and cools radia- 
tively, the resulting X-ray luminosity can be approximated 
as 
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where the term in bracketts is the stellar mass loss rate 
per blue luminosity (cf. Faber & Gallagher 1976; Sarazin 
1990). For NGC 4550, with cr* w 300 km s" 1 and L B = 
0.3 x 10 10 Lq, we get L x « lx 10 39 ergs s _1 . This is 
consistent with the existing upper limit for NGC 4550 of 
4 x 10 40 ergs s" 1 (Fabbiano et al. 1992). The time scale for 



accumulating an X-ray gas mass equivalent to the molec- 
ular gas mass is 



At 



Lb 



L^M g 



300 Myr . 



(4) 



[M*/Lb] is the stellar mass loss rate from evolved stars, 
assumed to be 1.5 x 10~ n M Q year -1 Lq 1 (cf. Faber k, 
Gallagher 1976; Sarazin 1990). 

Hence, a substantial hot X-ray emitting plasma can be 
established in the center of NGC 4550 on a relatively short 
time scale. This plasma can contribute to the heating of 
the dust component as well as contribute to its destruction 
through sputtering. 

Assuming a core radius for the X-ray gas of 1 kpc 
(12" at D = 16.8Mpc), the average X-ray intensity in 
the central region of NGC 4550 is I\ — 3 x 10 -6 ergs s _1 



resulting in a heating rate per dust grain of 
1 x 10~ 14 ergs s _1 . This is three orders of 



cm -2 ster -1 , 
47r lx 7r a 2 R: 

magnitude smaller than the heating rate from the ambient 
stellar light (Sect. 3.3). The X-ray photon contribution to 
the heating of the dust can thus be neglected. 

The hot electrons in the core of the X-ray emitting gas 
will collide with dust grains. The flux of hot electrons is 
F e = ^(SfcTe/Trme) 1 / 2 = 6.2 x 10 8 n e (T e /10 6 K) 1 / 2 cmT 2 
s _1 (cf. de Jong et al. 1990). Most of the kinetic energy 
of the hot electrons will be deposited in the grains upon 
collision, resulting in a heating rate (3/2) F e k T e ira 2 m 
4.0 x 10~ n n e (T e /10 e K) ergss -1 . Assuming an electron 
density of 10~ 2 cm -3 and the above derived gas temper- 
ature of 2 x 10 6 K, the heating rate per dust grain due 
to hot electrons becomes 1 x 10~ 12 ergss -1 . This is much 
larger than the heating rate due to X-ray photons but still 
ten times less than the heating rate due to the ambient 
stellar radiation field. 

The hot plasma will destroy the dust grains through 
'sputtering' (cf. Draine & Salpeter 1979). For diffuse gas 
(i.e. Ay < 1, the time scale is very short: At w 2 x 
10 5 (rig/cur 3 )- 1 (a/0.1/xm)yr (Draine & Salpter 1979). 
With n e w 10~ 2 cm -3 and a typical grain size of 0.1/im, 
At is only 20 Myr. After ~100 Myr the only surviving gas 
clouds are those which are optically opaque with Ay > 3 
(cf. de Jong et al. 1990). The lifetime of these opaque 
clouds depends strongly on geometry and physical condi- 
tions. De Jong et al. (1990) estimated typical parameters 
for this dense gas component in NGC 4696, the central el- 
liptical galaxy of the Centaurus cluster, and found values 
similar to the cores of molecular clouds in our Galaxy. Star 
formation could be sustained in these environments. 

To summarize: Mass lost from evolved stars will be 
heated to a temperature of ~ 2 x 10 6 K on a relatively 
short time scale due to the high (one dimensional) stellar 
velocity dispersion in NGC 4550. The X-ray luminosity is 
consistent with the upper limit from Einstein data. The 
X-ray gas will not contribute significantly to the heating 
of the dust grains but will destroy any diffuse component 
through sputtering on a time scale of 20-100 Myr. The 
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only remaining dusty molecular gas component is associ- 
ated with dense clouds of high optical opacity, possibly 
sustaining star formation. 

3.5. Star formation 

Assuming for the moment that the FIR luminosity is pow- 
ered entirely by young and massive stars and that the ini- 
tial mass function (IMF) is of Salpeter type with cut-offs 
at 0.1 and 100 M Q (cf. Thronson & Telesco 1986), the 
star formation rate (SFR) becomes 0.065 M Q yr _1 . The 
ratio of Lfir/A/h 2 is usually taken as a measure of the 
star formation efficiency. In NGC 4550 this ratio is 8 ± 2 
Lq/M©, which is similar to the values found for normal 
spiral galaxies. The gas consumption time scale is defined 
as the ratio Af H2 /SFR, which for NGC 4550 is 2 x 10 8 
years. 

Is star formation really going on in the center of 
NGC 4550? Kennicutt (1989) suggested that significant 
star formation only occurs only when the gas surface den- 
sity S gas exceeds a critical value Sg"s- The latter value 
depends on the galaxy rotation and on the gas velocity 
dispersion. The rotation curve presented by Rubin et al. 
(1992, 1997) shows that the ionized gas rotates approx- 
imately like a solid body in the inner 7", with a maxi- 
mum rotational velocity of 160 kms -1 . Within the region 
of solid body rotation we can express the critical gas sur- 
face density as 

S™* = 0.138afic M pc -2 , (5) 

where ft is the angular velocity in units of kms -1 kpc -1 , 
c is the gas velocity dispersion in units of kms" 1 and 
a is a dimensionless constant of order unity. For solid 
body rotation, this relation is independent of galactocen- 
tric distance as long as the velocity dispersion remains 
constant. We do not know the extent of the molecular gas 
disk nor its velocity dispersion, but reasonable assump- 
tions are that it coexists with the dust (Fig. |3j), having 
a radial extent of ~7", and that it has a velocity disper- 
sion ~10kms . Furthermore, assuming that the molecu- 
lar gas corotates with the ionized gas, the critical gas sur- 
face density is Egj* « 260(c/10kms" 1 ) M pc" 2 . The ob- 
served gas surface density, assuming a homogeneous disk 
with a radial extent of 7" (0.58 kpc) is £ gas « 17 M Q pc -2 , 
where we have included He at primordial abundances. The 
ratio Egas/E™* = 0.07(c/10kms -1 ) -1 . Since we found 
the molecular gas to be distributed non-axisymmetrically 
in the same sense as the dust seen in Fig. [| the ac- 
tual gas surface gas density is likely a factor of two 
larger. Nevertheless, even if the velocity dispersion is only 
5 kms -1 , the gas surface density is only ~ 0.25E™*. This 
means that the molecular gas is likely to be stable against 
gravitational collapse on a global scale. Smaller scales with 
potentially higher gas surface densities will be discussed 
in the following section. 



4. Discussion 

Although NGC 4550 is exceptional in that it contains two 
similar stellar disks counterrotating with respect to each 
other, there are few if any morphological and photomet- 
rical properties that distinguish it from other early-type 
galaxies. Surface photometry of 14 E and SO galaxies in the 
Virgo cluster, obtained with the Hubble Space Telescope 
(HST) and including NGC 4550 ( Jaffe et al. 1994) reveals 
nuclear disks in all of the early type galaxies later than 
E4. A nuclear dust component was also seen in most 
of the early type galaxies. Jaffe et al. (1994) did not 
find any property distinguishing NGC 4550 from the other 
13 galaxies, in spite of the counterrotating stellar disks. 
Furthermore, the FIR properties of the early type galax- 
ies in the Jaffe et al. (1994) sample are similar to those of 
NGC 4550. Hence, the presence of ~ 10 7 M Q of molecular 
gas in the center of NGC 4550 is not exceptional for its 
morphology (see also van Dokkum & Franx 1995). The 
only peculiarity which may point to the presence of the 
counterrotating disks is the non-axisymmetric distribution 
of the molecular gas and dust in the center of NGC 4550. 

Several numerical/theoretical studies of the dynam- 
ics of disks containing counterrotating stars and/or gas 
has shown that the fastest growing unstable mode is an 
to — 1 instability (Merrifield & Kuijken 1994; Lovelace 
et al. 1997; Thakar et al. 1996; Thakar & Ryden 1998; 
Garci'a-Burillo et al. 2000). This applies to the cases where 
the mass fractions of the pro- and retrograde mass compo- 
nents are of comparable size. The result is a tightly wound 
one-armed spiral structure within a few hundred pc from 
the nucleus. Initially the to = 1 instability can also be 
manifest through a lopsidedness of the mass distribution. 
This latter feature is, however, relatively short-lived, 0- 
500 Myr. On time scales longer than 1 Gyr an to = 2 
instability dominates (cf. Garcfa-Burillo et al. 2000). This 
takes the form of a bar or oval configuration. 

An extensive study of instabilities in counterrotating 
self-gravitating stellar disks was presented by Sellwood & 
Merritt (1994). They found that for kinematically cool sys- 
tems, lopsidedness is the only instability which is unique 
for counterrotating stellar systems (this applies for iso- 
lated galaxies). Contrary to the studies mentioned above, 
Sellwood & Merritt found the lopsidedness to persist for 
long periods and suggested that this feature could be used 
to identify galaxies with counterrotating stellar subsys- 
tems. 

The dust distribution in NGC 4550 is reminiscent of a 
tightly wound spiral arm structure, extending to a galac- 
tocentric distance of ~600pc, but with a pronounced lop- 
sidedness. As described in Sect.||, the CO emission is con- 
sistent with an asymmetric distribution of the molecular 
gas, in the same manner as the distribution of the dust 
component. The kinematics of the CO emission suggests 
that molecular gas is present on both sides of the nucleus, 
but that the northern side is dominant. In view of the 
above mentioned numerical simulations, the most likely 
explanation for the lopsided gas and dust distribution in 
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Fig. 3. HST WFPC2 images of the central region of NGC 4550. Only the PC chip is shown, with a pixel size of 0.046". 
Three images are displayed, centered on the nucleus: F555W (V-band), F814W (I-band) and the V-I image (inverted). 
North is up and east to the left. 



NGC 4550 is an instability caused by the counterrotating 
disks. As pointed out by Sellwood & Merritt (1994), the 
lopsided instability is strongest for the coolest systems, 
i.e. the gas. It remains to be shown whether these features 
can be longlived. If not, the lopsidedness implies a recent 
accretion event (<0.5 Gyr). 

In Sect. 3.4 we argued for the existence of a hot X- 
ray emitting plasma in the central region of NGC 4550. 
Although this X-ray component was found to give a 
neglible contribution to the heating rate of the dust grains, 
it limits the lifetime of any diffuse dust component to 
less than ^100 Myrs. A dense gas component can sur- 
vive for longer time scales. The physical characteristics 
needed for molecular gas surviving in a hot X-ray emit- 
ting plasma implies small and dense molecular clouds, sim- 
ilar to molecular cloud cores in in our own Galaxy. These 
cores could very well sustain star formation despite that 
the global surface density of molecular gas in NGC 4550 
was shown to be too low to be unstable to gravitational 
collapse (Sect |3~5| ). This star formation activity would con- 
tribute the extra heating suggested to explain the high 
dust temperature of 39 ± 7 K (see Sect. j3~3|). 

The only other known system with a significant coun- 
terrotating stellar disk and molecular gas is NGC 3593. 
In contrast to NGC 4550, this edge-on SO/a galaxy con- 
tains two stellar disks with different central surface bright- 
ness and different scale heights (Bertola et al. 1996). The 
smaller disk, which makes up ~18% of the total disk stellar 
mass, also contains significant amounts of ionized, atomic 
as well as molecular gas (Wiklind & Henkel 1992; Corsini 
et al. 1998; Garcfa-Burillo et al. 2000). The atomic and 
molecular gas masses are 1.8 x 10 s and 4.5 x 10 s M Q , re- 
spectively (Krumm & Salpeter 1979; Wiklind & Henkel 
1992). Including the counterrotating gas components, the 
mass ratio of the two disks is ~0.3, with the gas making 
up ^40% of the total mass of the counterrotating disk 
(Mstars + Mgas)- This is in stark contrast to NGC 4550, 
where the mass ratio of the two disks is ~1 with the gas 
making up ~0.1% of the mass of one of the disks. 

NGC 3593 shows several signs of vigorous star for- 
mation activity (cf. Henkel & Wiklind 1997). It has a 
relatively large FIR luminosity suggesting an SFR of 
l.eMoyr" 1 (Wiklind & Henkel 1992), 25 times higher 
than NGC 4550. The star formation activity is confined to 
a central ring with a radius of ~300pc. The ring is seen in 
both ionized and molecular gas (Wiklind & Henkel 1992; 
Corsini et al. 1998; Garcfa-Burillo et al. 2000). The atomic 
gas has only been observed with a single dish telescope but 
appears to be significantly more extended than both the 
ionized and molecular gas component. It is thus possible 
that in NGC 3593 we are witnessing the transformation 
of accreted gas, with an angular momentum opposite to 
that of the main stellar disk, into stars. Once the gas is 



consumed, NGC 3593 will look like an SO galaxy with a 
counterrotating stellar population comprising ~25% of the 
total stellar disk mass. Although less than in NGC 4550, 
it is still significant in comparison to the limits obtained 
for counterrotating stellar components in the sample of 28 
SO galaxies by Kuijken et al. (1996). 

In the case of NGC 4550 there is only weak and indi- 
rect evidence for on-going star formation. Regardless of 
whether star formation is taking place or not, the age of 
the observed molecular gas and dust must be low. Since 
mass loss from evolved stars would provide cool gas asso- 
ciated with both stellar disks that are old and approx- 
imately coeval (see Rix et al. 1992), the gas (and the 
dust) must be the result of a relatively recent (<Cl Gyr) 
accretion event not related to the main build-up of the 
two counterrotating stellar populations. Searching for po- 
tential sources for a mass transfer, the nearby elliptical 
galaxy NGC4551 (V sya ~ 1170 kms" 1 versus 380 kms" 1 
for NGC 4550) is a highly unlikely target; the accretion of 
a dwarf galaxy is therefore a more reasonable scenario. 

Acknowledgements. TW acknowledges support from NFR (the 
Swedish Natural Science Research Council) for a research 
grant. 



References 

Bertola, F., Buson, L. M., Zeilinger, W. W. 1992, ApJ 401, L79 
Bertola, F., Cinzano P., Corsini E. M., et al., 1996, ApJ 458, 
L67 

Corsini, E. M., Pizzella, A., Funes, J. C, Vega Beltran, J. C, 

Bertola, F. 1998, A&A 337, 80 
de Jong, T., N0rgaard-Nielsen, H. U., J0rgensen, H. E., 

Hansen, L. 1990, A&A 232, 317 
Draine, B. T., Salpeter, E. E. 1979, ApJ 231, 77 
DuPrie, K., Schneider, S. E. 1996, AJ 112, 937 
Faber, S. M., Gallagher, J. S. 1976, ApJ 204, 365 
Fabbiano, G., Kim, D.-W., Trinchieri, G. 1992, ApJS 80, 531 
Faundez, S., Kenney, J., Murphy, E. 2000, AAS 196, #29.15 
Garcfa-Burillo, S., Sempere, M. J., Combes, F., Hunt, L. K., 

Neri, R. 2000, A&A 363, 869 
Henkel, C, Wiklind, T. 1997, Sp. Sc. Rev. 81, 1 
Hildebrand, R. H. 1983, QJRAS 24, 267 

Ho, L. C, Rudnick, G., Rix, H.-W., et al. 2000, ApJ 541, 120 
Huchtmeier, W. K., Richter, O.-G., 1989, A General Catalog 

of H I Observations of Galaxies, Springer Verlag, New York 
Jaffe, W., Ford, H. C, O'Connell, R. W., van den Bosch, F. 

C, Ferrarese, L. 1994, A J 108, 1567 
Jones, T. W., Merrrill, K. M. 1976, ApJ 209, 509 
Jore, K. P., Broeils, A. H., Haynes, M. P. 1996, AJ 112, 438 
Kennicutt, R. C. 1989, ApJ 344, 685 
Krumm, N., Salpeter, E. E. 1979, ApJ 228, 64 
Kuijken, K., Fisher, D., Merrifield, M. R. 1996, MNRAS 283, 

543 

Lovelace, R. V. E., Jore, K. P., Haynes, M. P. 1997, ApJ 475, 
83 



8 



Wiklind & Henkel: Molecular gas and dust in NGC 4550 



Merrifield, M. R., Kuijken, K. 1994, ApJ 432, 575 

Pfenniger, D. 1999, IAU Symp. 186, in Galaxy Interactions at 

Low and High Redshift, eds. J.E. Barnes, D.B. Sanders, 

Kluwer Academic Publishers, Dordrecht, p. 157 
Prada, F., Gutierrez, C. M., Peletier, R. F., McKeith, C. D. 

1996, ApJ 463, L9 
Rix, H.-W., Franx, M., Fisher, D., Illingworth, G. 1992, ApJ 

400, L5 

Rubin, V. C, Graham, J. A., Kenney, J. D. P. 1992, ApJ 394, 
L9 

Rubin, V. C, Kenney, J. D. P., Young, J. S. 1997, AJ 113, 
1250 

Sandage, A., Visvanathan, N. 1978, ApJ 223, 707 

Sarazin, C. L. 1990, in The Interstellar Medium in Galaxies, 

eds. H. A. Thronson, J. M., Shull, Kluwer Academic 

Publishers, Dordrecht, p. 201 
Sellwood, J. A., Merritt, D. 1994, ApJ 425, 530 
Thakar, A. R., Ryden, B. S. 1996, ApJ 461, 55 
Thakar, A. R., Ryden, B. S., Jore, K. P., Broeils, A. H. 1997, 

ApJ 479, 702 
Thakar, A. R., Ryden, B. S. 1998, ApJ 506, 93 
Thronson, H. A., Telesco, C. M. 1986, ApJ 311, 98 
Tully, R. B. 1988, Nearby Galaxies Catalog, Cambridge 

University Press, Cambridge 
van Dokkum, P. G., Franx, M. 1995, AJ 110, 2027 
Wiklind, T., Henkel, C. 1992, A&A 257, 437 



This figure "ngc4550-fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0106286vl 



